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Abstract 

We perform molecular dynamics simulations to investigate the reduction of the thermal conduc- 
tivity by kinks in silicon nanowires. The reduction percentage can be as high as 70% at room 
temperature. The temperature dependence of the reduction is also calculated. By calculating 
phonon polarization vectors, two mechanisms are found to be responsible for the reduced thermal 
conductivity: (1) the interchanging effect between the longitudinal and transverse phonon modes 
and (2) the pinching effect, i.e a new type of localization, for the twisting and transverse phonon 
modes in the kinked silicon nanowires. Our work demonstrates that the phonon interchanging 
and pinching effects, induced by kinking, are brand new and effective ways in modulating heat 
transfer in nanowires, which enables the kinked silicon nanowires to be a promising candidate for 
thermoelectric materials. 

PACS numbers: 62.23.Hj, 65.80.-g, 63.22.-m, 63.22.Gh 
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I. INTRODUCTION 



In 2009, a new type of building block - kinked silicon nanowire (KSiNW) - was synthe- 
sized by Tian et al. in Leiber's group at Harvard University- Particularly, they were able 
to manipulate the arm length of the kink by controlling the growing time. Since then, great 
experimental efforts have been devoted to investigating the growth mechanism of the kinks 
in silicon nanowires (SiNWs) Besides these experimental works, Schwarz and Tersoff have 
proposed a theoretical model to interpret the growth mechanism of KSiNWs. In their model, 
the formation of kinks is due to the interplay of three basic processes: facet growth, droplet 
statics, and the introduction of new facets.-^ Stimulated by KSiNWs, several groups have 
synthesized kinks in other nanowires, such as In 2 03 multi-kinked nanowires,— kinked ger- 
manium nanowires,— germaniumsilicon axial heterostructure with kinks,— and kinked ZnO 
nanowires.— While existing works for KSiNW mainly concentrate on its growth mechanism, 
a new experiment from Leiber's group has demonstrated the application of KSiNWs as na- 
noelectronic bioprobes.— However, the study of thermal properties like thermal conductivity 
in KSiNWs is quite limited, although this is important for its application in thermoelectrics. 

Recently, the thermal transport of SiNWs has attracted broad interest, because phys- 
ical understanding of the thermal properties in nano-materials remains unclear.— ~— Volz 
and Chen predicted the ultra-low thermal conductivity of SiNWs, which was about two 
orders of magnitude smaller than the thermal conductivity of bulk Si crystals obtained from 
molecular dynamics (MD) simulations.— A few years later, Li et al. showed experimentally 
the ultra-low thermal conductivity of single SiNW.— In 1998, the phonon confinement was 
firstly proposed as a mechanism for the thermal conductivity reduction in low- dimensional 
structures by Balandin and Wang,— and this idea was worked out specifically to SiNWs 
and Si-Ge NWs in the following years.— Ponomareva et.al also observed this phonon 
confinement effect on the thermal conductivity in SiNWs in 2007.— Yang et.al found very 
low thermal conductivity in the isotope-doped SiNWs.— By introducing surface roughness, 
it was found that the thermal conductivity in nanowires can be greatly reduced from both 
theoretical^ 1 ^ and experimental examinations^. For pristine SiNWs, the thermal conduc- 
tivity is found to be strongly anisotropic— Due to the effects of boundary scattering and 
confinement on lattice thermal conductivity, there is a significant progress on the applica- 
tions of SiNWs. For example, SiNWs can be used in thermoelectric systems where ultra-low 
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thermal conductivity materials are needed. 

In this letter, we perform non-equilibrium MD simulations to study the thermal conduc- 
tivity in KSiNWs and compare it to that of straight silicon nanowires (SSiNWs). We find 
that the thermal conductivity of SSiNWs can be reduced up to 70% by kinks at room tem- 
perature and the reduction percentage depends on the temperature. After comparatively 
studying the phonon modes in the SSiNW and the KSiNW, we show that the reduction in 
the thermal conductivity comes from the kink-induced interchanging and pinching effects 
on some phonon modes in the KSiNW. 

II. STRUCTURE AND CALCULATION DETAILS 

Fig. [1] displays the configuration of three KSiNWs of 11.0 A in diameter and 190.0 A in 
length along the axial direction. The growth direction changes from one < 110 > to another 
< 110 > direction at the kink. The angle between these two < 110 > directions is 90°. 
Fixed boundary condition is applied in the growth direction and free boundary condition 
is applied in the lateral direction. The force between silicon atoms in the Newton equation 
is described by the commonly used Tersoff bond order potentials.— The Newton equations 
are solved by the velocity Verlet algorithm with an integration time step of 1.0 fs. A total 
simulation time is typically around 2.0 ns and will be sufficiently extended to guarantee the 
achievement of steady state whenever necessary. 

The two regions close to the left/right ends are maintained at high/low temperatures 
Tl/r = (1± q)T by Nose-Hoover thermostat.— 1 ^ T is the average temperature, a = 0.1 
is adopted in our calculation. The heat energy is pumped into the system through the 
left temperature-controlled region with thermal current Jl and will flow out of the right 
temperature-controlled region with current Jr. At steady state, we have Jl = —Jr as 
required by the energy conservation law. Using this relation, the thermal current across the 
system can be obtained by J = (Jl — Jr)/2, with dJ = Jl + Jr as an estimated error.— 
The thermal conductivity is obtained from the Fourier law. There is no quantum correction 
here, as all simulations are performed above room temperature. The quantum correction 
will become more important for lower temperature. For the lowest temperature in present 
work, i.e 300 K, the quantum correction can be estimated to be around 3.0%, following the 
approach suggested by Maiti, Mahan, and Pantelides.— 
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III. RESULTS AND DISCUSSION 



Fig. [2] shows the relative thermal conductivity versus the number of kinks of KSiNWs. 
The relative thermal conductivity of 1 corresponds to the value of thermal conductivity for 
the SSiNW, which are 7.9 and 5.8 Wm" 1 K~ 1 at 300 and 600 K, respectively. All systems 
have the same diameter (D = 11.0 A) and the same length {L = 190.0 A) along the axial 
direction. We have observed similar reduction in KSiNWs with other diameter and length. 
As expected, the reduction of the thermal conductivity is obvious in the kinked structures. 
The relative thermal conductivity decreases to 44% and 31% as the kink number increases 
from to 11 at temperature 300K and 600 K, respectively. However, the relative thermal 
conductivity increases when keep increasing the kink number. In our simulation structures, 
the arm length Z arm (see Fig. [T|) is shortened as the kink number increases, which may be 
convenient for phonon transport. Similarly, there is also a minimum thermal conductivity 
when increasing the interface number in superlattice structures.— For ballistic heat 
transport, the minimum thermal conductivity results from the competition between the 
number and the localization properties of the confined phonon modes.— For diffusive heat 
transport, the minimum thermal conductivity is attributed to the interplay between the 
wave and the particle nature of the phonon quasi-particles.— 

As shown in Fig. |2l the reduction of the thermal conductivity also depends on the temper- 
ature. To investigate the temperature effect further, we show the temperature dependence of 
the thermal conductivity of both SSiNWs and KSiNWs in Fig. [3] (a). The thermal conduc- 
tivity in SSiNW decreases as the temperature increases because phonon-phonon scattering 
dominate the heat transport at high temperature. However, the thermal conductivity in 
KSiNW is insensitive to the temperature. Fig. [3] (b) shows the thermal conductivity ratio 
of KSiNW to SSiNW based on the data in Fig. [3] (a). The relative thermal conductivity 
increases monotonically with increasing temperature, because the thermal conductivity in 
SSiNWs is largely deduced at high temperature and the thermal conductivity in KSiNWs 
is temperature insensitive. There are less phonon-phonon scatterings and a higher thermal 
conductivity in SSiNWs at low temperature, where the phonon pinching is more effective in 
confining the heat transfer. 

To reveal the microscopic origin for the reduction of the thermal conductivity in the 
KSiNW, we comparatively study the phonon modes in the SSiNW and KSiNW, which are 



4 



the only energy carrier in the non-metal solids. In Fig. HJ we show the dispersion of all four 
acoustic branches of SSiNW—, which are one longitudinal acoustic (LA), one twisting (TW), 
and two transverse acoustic (TA) modes. The full phonon dispersion of SSiNW is shown in 
the right bottom inset of Fig. HI 

To show the pinching effect in modes, we calculate the eigen modes of the four acoustic 
branches in the SSiNW with length L = 190 A and a KSiNW with one kink and 95 A in 
arm length. The diameter is kept the same for all structures. Both ends of the nanowire are 
fixed in this calculation, which is consistent with the fixed boundary condition applied in 
the above thermal transport simulation. The phonon modes (frequencies and polarization 
vectors) are obtained by solving the eigenvalue problem of the dynamical matrix derived 
from the Tersoff potential.— The dynamical matrix is obtained by = cPV/dxidXj, where 
V is the Tersoff potential and X{ is the position of the i-th degree of freedom. This formula 
is realized numerically by calculating the energy change after a small displacement of the 
i-th and j-th degrees of freedom. The obtained phonon modes from this calculation for the 
finite system are related to the traveling phonons in Fig. H]for a bulk structure (with periodic 
boundary condition applying on a small unit cell). The former is essentially a special solution 
of the latter under particular boundary conditions (fix boundary conditions here), so the 
phonon modes in the finite system can be decomposed into some traveling waves which will 
change into the corresponding reflection and tunneling waves at the kinks. As a result, the 
phonon modes in the finite system provide intrinsic information for the thermal transport. 

In the following, the eigen modes in the SSiNW are compared to those in the KSiNW. 
Four modes with the lowest frequency in each acoustic branch in SSiNWs are selected and 
presented in Fig. |5j The polarization vectors of the eigen-modes are represented by arrows 
whose length corresponds to the amplitudes. Red (blue) arrows indicate a mode with good 
(poor) heat transport capability; i.e the phonon is extended (pinched). To analyze the 
correspondence between phonon modes from both structures in Fig. |5l we have examined 
the polarization vectors of the first thirty lowest-frequency phonon modes in the SSiNW and 
the KSiNW. 

Fig. |5] (a) shows the polarization vectors of a LA mode (along x direction) with 
u> = 7.4 cm -1 in the SSiNW. In the KSiNW, the polarization vector of a LA mode with 
u = 7.7 cm -1 is not confined in x direction. When the phonon transfers along the KSiNW 
from one end to the other, the LA mode and TA y mode interchange into each other, which 
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leads to the lose of the coherence of the phonon mode and eventually reduces its heat trans- 
port capability. This deduction effect is also observed in the folded graphene nanoribbon 
structure.— The frequency of the LA mode in the KSiNW is only slightly higher than that 
in the SSiNW due to the kink effect. We note that there is a mirror reflection symmetry 
locates at the middle of the kin angle in the KSiNW. Due to this mirror symmetry, there is 
another LA mode in the KSiNW with closer frequency (8.0 cm -1 ). 

In Fig. [5] (b), we compare a TA Z mode in the SSiNW and the KSiNW, whose polarization 
vectors are all along the ^-direction. The TA Z mode is almost not affected by the kink 
interface, owing to its special vibrating morphology which is perpendicular to the kink (xy 
plane). There is barely a difference in the frequency and vibration between the TA Z modes 
in SSiNW and that in KSiNW. Hence, the TA Z modes are good at heat delivery in both 
SSiNWs and KSiNWs. 

In the following, the phonon pinching is shown obviously in TA y and TW modes of 
KSiNWs. As shown in Fig. [5] (c), there is a critical difference between the TA y mode in 
SSiNW and that in KSiNW. The TA y mode with u = 0.7 cm" 1 in the SSiNW is a quite 
extended mode. While in the KSiNW, it is interesting that vibrating amplitudes of the TA y 
mode (u = 2.8 cm" 1 ) are nearly zero at the kink. In other words, the polarization vector of 
this phonon mode is pinched-off by the kink. The polarization vectors of the TA y modes are 
symmetric about the mirror plane at the kink. There is another TA y with a closer frequency 
(2.0 cm -1 ) in the KSiNW, whose polarization vectors are anti-symmetric about the mirror 
plane, and the phonon pinching effect is also found there. 

The pinching-off picture can be validated by the relationship between frequencies of 
corresponding modes in SSiNWs and KSiNWs. The frequency of the mode in the KSiNW 
(2.8 cm -1 ) is four times of the frequency in the SSiNW (0.7 cm -1 ). The arm length of the 
KSiNW is half of the length of the SSiNW. As shown in Fig. HI the frequency of the TA y 
mode is proportional to the square of the wave vector, oj = (3k 2 = (3(2tt/L) 2 , due to the 
flexure nature of the quasi-one-dimensional configuration of the nanowire,— /3 is a coefficient 
determined by the Tersoff potential. If a mode in SSiNW is pinched-off from the middle, 
the transverse mode confined in the left or right arm is degenerated, whose frequencies are 
/S(7r/L) 2 , i.e 4oj. This is the origin of the pinched-off transverse mode in the KSiNW. Similar 
to the TA y modes in the KSiNW, the pinching effect is observed in the TW modes shown 
in Fig. (d). Since the frequency for the TW mode is inversely proportional to the length 
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of the nanowire (Fig. H]), the frequency of the TW modes in the KSiNW is almost twice the 
frequency in the SSiNW. 

As a result of the pinching effect, the heat energy carried by the TA y and the TW modes 
is difficult to be transported from the left arm to the right arm of the KSiNW. The pinched 
phonon mode is a new type of localization mode. This is the main reason for the reduction 
of the thermal conductivity by kinks in the KSiNW. The dependence of the reduction in 
the thermal conductivity of a KSiNW with fixed axial length (L) on kinks shown in Fig. [2] 
can be understood as follows: on the one hand, there are more pinching effects in KSiNW 
with more kinks, so the thermal conductivity becomes smaller with increasing kink number; 
on the other hand, the arm length is shorter for KSiNW with more kinks, and the pinching 
effect is weaker in KSiNW of shorter arm length as shown in Fig. O The arm length for the 
KSiNW studied in Fig. © is I 

arm — 47.5 A. We make a comparison between the frequencies 
in the SSiNW (w s ) and the KSiNW (co K ). It is obvious that w£A < 4wJ A and w£ w < 2wJ w , 
for the TA y and the TW modes, respectively. This provides a direct evidence for the weaker 
pinching effect in KSiNW with shorter arm length, because it is expected that ui^ A ps 4cjg A 
and Co>k W ~ w in case of a strong pinching effect such as that in Fig. Related to 
this weaker pinching effect, the vibrating amplitude of the atoms in the kink region starts 
to arise, which induces a tunneling phenomenon of this mode at the kink. This tunneling 
effect facilitates the heat delivery across the kink, thus enhancing the thermal conductivity 
of the KSiNW. The interplay between the pinching enhancement effect (due to increasing 
kink number) and the tunneling enhancement effect (due to the shortening of arm length) 
yields the minimum thermal conductivity in Fig. [21 

We would like to stress that lower thermal conductivity in nanowires is essential for 
their application in thermoelectric materials. There are some other methods to reduce the 
thermal conductivity, besides the phonon pinching effect discussed in present work. For 
example, Nika et.al demonstrated the importance of the phonon filtering effect on the ther- 
mal conductivity in one-dimensional superlattice nanowires.— They found very low thermal 
conductivity, because significant number of phonon modes are filtered out from thermal 
transport. The thermal conductivity can be further manipulated by combing the phonon 
pinching and other effects (eg. phonon filtering) in nanowires. 
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IV. CONCLUSTION 



We have performed MD simulation to study the heat transport in the KSiNW. We ob- 
serve that the thermal conductivity in the KSiNW is insensitive to the temperature, since 
the heat transport in KSiNW is mainly limited by the phonon scattering at the kink in- 
terface. We found that the thermal conductivity can be reduced by kinks in the KSiNW 
for 70% compared to that of the SSiNW at room temperature and the reduction can be 
further increased at lower temperature. By comparatively investigating the phonon modes 
in SSiNW and KSiNW, we disclose the underlying mechanism for the reduction of the ther- 
mal conductivity to be the kink-induced interchanging effect between LA and TA y modes 
and the pinching effect on TA y and the TW modes in the KSiNW. Both interchanging and 
pinching effects greatly suppress the transport ability of these phonon modes, leading to the 
reduction of the thermal conductivity. 
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FIG. 1: Configuration for KSiNWs with diameter D = 11.0 A and axial length L = 190.0 A. I, 
is the arm length. The number of kinks are 5, 7, and 11 in structures from top to bottom. 
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FIG. 2: Relative thermal conductivity versus number of kinks in KSiNWs of 11.0 A in diameter 
and 190.0 A in length along the axial direction. The value of thermal conductivity for the SSiNW 
are 7.67 and 5.78 Wm~ 1 K~ 1 at 300 and 600 K, respectively. 
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FIG. 3: Thermal conductivity versus MD simulation temperature, (a) Thermal conductivity in 
SSiNW (dots) and KSiNW with 11 kinks (pentagons), (b) The ratio of thermal conductivity in 
KSiNW to that in SSiNW. The lines are guides to the eye. 




FIG. 4: Phonon dispersion of the four low-frequency acoustic modes: LA, TW, and two TA modes 
in SSiNW of 11.0 A in diameter and period length A = 33.0 A. Left top inset gives the cross-section 
of the nanowire. Right bottom inset shows the full phonon dispersion in the first Broullin zone. 
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FIG. 5: (Color online) Comparison between the four low-frequency acoustic modes in the SSiNW 
with length L = 190 A and the KSiNW with arm length Z arm = 95 A. Arrows in the figure represent 
the polarization vector of the phonon mode. Red (blue) color indicates good (poor) heat transport 
capability for the mode, (a) LA mode, (b) TA mode vibrating in z direction, (c) TA mode 
vibrating in y direction, (d) TW mode. Numbers in the figure are the corresponding frequency in 
unit of cm -1 . Note the pinching effect for the two modes in the KSiNW in (c) and (d). 
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FIG. 6: (Color online) Phonon pinching effect in the SSiNW with length L = 95 A and the 
KSiNW with arm length Z arm = 47.5 A for (a) TA mode vibrating in y direction, and (b) TW 
mode. Pinching effect here is weaker than the two corresponding modes in the KSiNW of longer 
arm length shown in Fig. [5] (c) and (d). 
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